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Abstract 

We explore the device potential of tunable-gap bilayer graphene FET exploiting the possibility of opening a 
bandgap in bilayer graphene by applying a vertical electric field via independent gate operation. We evaluate device 
behavior using atomistic simulations based on the self-consistent solution of the Poisson and Schrodinger equations 
within the NEGF formalism. We show that the concept works, but bandgap opening is not strong enough to suppress 
band-to-band tunneling in order to obtain a sufficiently large I on /I s ratio for CMOS device operation. 
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I. Introduction 

Graphene has entered the nanoelectronics scenario only recently and it is intensely investigated as one of the 
most promising candidates to replace silicon as a channel material in nanoscale transistors. Even though graphene 
is a gapless material, a significant energy gap can be opened in different ways, like for example by "rolling" it in 
carbon nanotubes [1] or by defining narrow graphene stripes through electron beam lithography [2] or chemical 
reaction [3]. However, several unsolved technological problems arise: state of the art technology cannot indeed 
conveniently control the chirality and the nanotube position, as well as define graphene nanoribbons with widths 
close to 1-2 nm with atomically flat edges, in order to obtain a reasonable semiconducting energy gap. 

Recently, theoretical models [4] and experiments [5], [6] have shown that bilayer graphene has an energy gap 
controllable by a vertical electric field. One could exploit this property to use bilayer graphene as a channel material 
for FETs, defining an energy gap only when really needed, i.e. when the FET must be in the off state. This capability 
could open the possibility of patterning a bilayer graphene sheet with lithographic techniques and non prohibitive 
feature sizes. 

In this letter, we want to assess the possibility of fabricating a tunable-gap Bilayer Graphene (BG) FET with 
independent gate operation, based on atomistic numerical simulations. We assume that all the technological chal- 
lenges associated to the reliable fabrication of bilayer graphene FETs can be solved, and evaluate the potential 
performance of near ideal device structures. In our opinion, this is one of the most powerful uses of computer 
simulations in guiding and orienting nanoelectronics research. 

A semi-analytical model based on the effective mass approximation has been presented in order to compare 
bilayer against monolayer graphene transistors [7]. However, in that work a constant induced energy bandgap is 
assumed for bilayer graphene, missing the most important material property of the material. 

For a physically sound analysis of BG-FETs, we have developed a code based on the Non-Equilibrium Green's 
Function formalism (NEGF), with a tight-binding Hamiltonian on a p z orbital basis set in the real space, which 
has been included in our in-house device simulator NANOTCAD ViDES [8]. 

We will show that BG-FETs miss the ITRS requirements [9] for the I on /I s ratio (> 10 4 ) by a large amount, 
since the induced gap is not sufficient to suppress band-to-band tunneling currents. 

II. Physical model and Results 

Our approach is based on the self-consistent solution of the Poisson and Schrodinger equations within the NEGF 
formalism and the ballistic transport assumption, as described in [10]. 

The considered bilayer graphene Hamiltonian is composed by the two single layer graphene Hamiltonians, 
coupled by the t p =0.35 eV hopping parameters in correspondence of overlaying atoms along the z-direction [4]: 
the elementary cell is depicted in Fig. QJ. Semi-infinite contacts have been modeled along the y direction by means 
of self-energies, while periodic boundary conditions have been imposed in the x direction, with period equal to 
\fia cc , where a cc = 0.144 nm is the carbon-carbon bonding distance. 



October 1, 2008 



DRAFT 



4 



The simulated device is a double-gate BG-FET, whose structure is shown in Fig. QJ). We assume metal gates, and 
a 1.5 nm layer of SiC>2 as gate dielectric. We also assume an air gap of 0.5 nm between the dielectric interface and 
the position of carbon sites, as evaluated in [3], The channel is 15 nm long, and the inter-layer distance is 0.35 nm. 
The source and drain extensions are 10 nm long, and are doped with an equivalent molar fraction of fully ionized 
donors f d = 5 x 10~ 3 . 

Particular attention has to be posed in the computation of the mobile charge. Considering for the sake of simplicity 
the equilibrium case, the mobile charge on a given site on layer j (j = 1, 2) of the bilayer graphene sheet reads 

r+oc 

Pj = -2q / dE LDOSj(E)f(E) 

J Ei 
r+oo 

+2q dE LDOSj(E){l- f(E)} (1) 

JEi 

where q is the elementary charge, Ei is the intrinsic (mid-gap) Fermi level, / is the Fermi-Dirac occupation factor, 

and LDOSj(E) is the local density of states on layer j. 

Let us stress the fact that, in order to avoid an unphysical antiscreening behavior, it is necessary to include the 

effect of dielectric polarization along the direction perpendicular to the graphene bilayer, due to the response of 

valence electrons to the applied field. We include this effect by introducing a polarization factor a defined as 

«n J^LDOS^E^dE 

<x(V) = —p , (2) 

[LDOSi(E, V) + LDOS 2 {E, V)] dE 

where V is the potential difference between the two graphene layers. By using a, we can write the total charge 
(mobile charge, valence band electrons and ions) on a given site of layers 1 and 2 as 

Pi = pi+q[l-2a(V)] 

p' 2 = p 2 -q[l-2a(V)}. (3) 

However, performing the integral from — oo to E^ in (O at each iteration step can be too computationally 
demanding. To this purpose, we compute the polarization factor as a function of V only once, for an infinitely 
long BG-FET (solid line in Fig. [it). As can be seen, a(V) can be reasonably approximated by a straight line 
a(V) = 0.5 + 3.282 x 10 -3 y. In our code, we have used an odd polynomial of 11th order obtained through least 
mean square fitting (symbols in Fig. [it). 

In Fig. [2] the transfer characteristics for Vds = 0.5 V are shown, for different gate configurations. In particular, 
in Fig. [2t the transfer characteristics as a function of the top gate voltage Vtop for two bottom gate voltages 
(Vbottom — V and Vbottom = —1-0 V) are shown. As can be seen, large currents can be obtained above 
threshold, while the BG-FET shows huge problems to switch off. Even after applying really negative gate voltages, 
the drain-to-source current decreases only by roughly six times. The ratio even worsens when Vbottom = —1-0 V, 
since in this case a larger amount of positive charge is injected in the valence band, due to band-to-band tunneling, 
which pins the potential along the channel, and further degrades the gate control over the barrier. 

Such results are close to those obtained in experiments on monolayer graphene devices [11], showing that the 
tunable gap is too small to effectively suppress band-to-band tunneling in the off state. 
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In order to increase the vertical electric field between the two layers of graphene, we have driven the device by 
imposing a fixed voltage between the two gates Vdig = Vt op — Vbottom- 

Even in this case, device transconductance (i.e. the slope of the transfer characteristic) almost remains the same, 
regardless of the applied Vdiff- This can be explained by the large accumulation of positive charge in the channel, 
caused by band-to-band tunneling in correspondence of the drain [10]: this charge screens the vertical electric field, 
limiting the size of the opened bandgap. 

In Fig [3]3, the transmission coefficient as a function of the transversal wave vector k x is shown for V top — —4.0 V, 
Vds = 0.5 V and Vbottom = V, i.e., close to minimum achievable current, and in correspondence of the band 
minimum. Different energy intervals are identified with letters and compared to a sketch of the band edge profile 
in Fig. [3^ as an aid to interpretation. As can be seen, energy intervals tagged as B, D and F correspond to the 
bandgaps in the channel, source and drain extensions, respectively: the gap is not uniform, it is maximum in the 
channel, where it is almost 100 meV. One can also see the bound states in the valence band in the energy interval 
C and the valence band states in E, which are populated by holes for negative voltages, and which eventually pin 
the potential of the two mono-layers of graphene, i.e. reducing the gate effectiveness in controlling the channel 
barrier and the bandgap. 

Band-to-band tunneling can in principle be reduced by engineering the electric field between source and channel. 
Here we make two attempts, by introducing a lateral spacer on the gate sides of 2.5 nm, and by introducing a 
linear doping profile degrading to zero in 2.5 nm under the spacer. Both options have a positive but limited effect, 
as shown in Fig. 3c, and do not solve the problem. 

III. Conclusions 

We have explored the possibility of realizing FETs by exploiting the tunable-gap property of bilayer graphene, 
through a code based on Tight-Binding NEGF device simulations. Aiming at technology foresight, we have 
considered the overly optimistic case of an ideal structure and of ballistic transport. 

Our computer simulations show that tuning the gap is not a very promising technique for achieving appropriate 
switching characteristics of BG-FETs, due to the too weak suppression of band-to-band tunneling. A more complete 
exploration of the design space, in terms of structure and bias, would be required to cast a definitive answer, that 
is beyond the scope of this letter. 
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Fig. 1. a) Elementary cell of the simulated bilayer graphene; b) transversal cross-section of the simulated graphene bilayer field effect transistor; 
c) tight-binding (solid line) and least mean square fitting analytical results (symbols) for the polarization factor a in eq. {2). 

Fig. 2. a) Transfer characteristics computed as a function of the top gate voltage, when the bottom gate voltage V^ ottom =0 V (solid line), and 
^6ottom = "l-0 V (dashed line); b) transfer characteristics computed as a function of Vtottom f°r Kj;//=1.0 V (solid line) and Vdiff=0.5 V 
(dashed line). 

Fig. 3. a) Sketch of the band edge profile; b) Transmission coefficient as a function of the energy and the transversal wave vector k x for 
^bottom— V, Vtop = —4.0 V and Vjjs = 0-5 V; b) transfer characteristics for a BG-FET with 2.5 nm spacer, and a BG-FET with doping 
profile linearly varying over a 25 nm long region from to fd- 
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FIG. 3 
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